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Measurement of the effective 



A precise determination of the effective B® — > K + K~ lifetime can be used to con- 
strain contributions from physics beyond the Standard Model in the B® meson 
system. Conventional approaches select B meson decay products that are signifi- 
cantly displaced from the B meson production vertex. As a consequence, B mesons 
with low decay times are suppressed, introducing a bias to the decay time spectrum 
which must be corrected. This analysis uses a technique that explicitly avoids a 
lifetime bias by using a neural network based trigger and event selection. Using 
1.0 fb _1 of data recorded by the LHCb experiment, the effective B® — > K + K~ 
lifetime is measured as 1.455 ± 0.046 (stat.) ± 0.006 (syst.) ps. 
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1 Introduction 



The study of charmless fe-hadron decays can be used to explore the phase structure of the 
CKM matrix and to search for indirect evidence of physics beyond the Standard Model 
(SM). A measurement of the effective lifetime of the B® — > K + K~ decay (charge conjugate 
modes are implied throughout) is of considerable interest as it is sensitive to new physical 
phenomena affecting the B® mixing phase and entering the decay at loop level [l}j4] . 

The B® — > K + K~ decay was first observed by the CDF collaboration (5] and the most 
precise measurement to date of the effective lifetime was made by the LHCb collaboration 



using data taken during 2010 



6j. A detailed theoretical description of the B® — > K + K 



decay can be found in Refs. (3]4j. When the initial flavour of the B® meson is unknown 
the decay time distribution can be written as 

T(t) ex (l-^Ar s )e~ r ^ + (l + ^ Ar Je- r ^. (1) 

The quantities T H and T L are the decay widths of the heavy and light B° s mass eigenstates 
and Ar s = IY — Th is the decay width difference. The parameter Aat s is defined as 
AAr a = — 2Re(A)/ (1 + |A| 2 ) where A = (q/p)(A/A), where the complex coefficients p and 
q define the mass eigenstates of the B° s -B° s system in terms of the flavour eigenstates (see 
e.g., Ref. jj]) and A (A) is the amplitude for a B Q S (B° s ) meson to decay to the K + K~ 
final state. 

If the decay time distribution given by Eq. [I] is fitted with a single exponential function 
the effective lifetime is given by [8] 



1 -vl 



1 + 2A A r a y s + yl 
1 + AAY a y s 



r B o(l + A Ars y s + 0(y 2 s )), (2) 



where t b o = 2/(T H + T L ) = Tj 1 and y s = Ar,/2r s . The K+R- final state is CP- 
even and so in the SM the decay is dominated by the light mass eigenstate such that 
*4Ar s = —0.972 ±0.012 IM9J and the effective lifetime thus is approximately equal to T L X . 



Adopting the approach of Ref. [3] and using the world averages of T s and Ar s 10 and 
the SM prediction of *4Ar s , the effective lifetime is predicted to be tkk = 1-40 ± 0.02 ps. 
However, the B®—> K + K~ decay is dominated by penguin diagrams and so is sensitive to 
physics beyond the SM entering at loop level, which may affect Aat s - The measurement 
is also sensitive to new physics contributions to the B® mixing phase which in turn affects 



Ar s 11 . Deviations from this prediction will therefore provide evidence of new physics. 

Conventional selections exploit the long lifetimes of 6-hadrons by requiring that their 
decay products are significantly displaced from the primary interaction point. However, 
this introduces a time-dependent acceptance of the selected 6-hadron candidates which 
needs to be taken into account in the analysis. This paper describes a technique based 
on neural networks which avoids such acceptance effects. Only properties independent 
of the decay time are used to discriminate between signal and background. To exploit 
the available information, including the correlations between variables, several neural 
networks are used in a dedicated trigger and event selection. 



1 



2 The LHCb experiment and simulation software 



The B° s — > K + K lifetime is measured using 1.0 fb 1 of pp collision data collected by 
the LHCb detector at a centre of mass energy of \/s = 7 TeV during 2011. The LHCb 



detector 12 is a single arm spectrometer with a pseudorapidity acceptance of 2 < rj < 5 
for charged particles. The detector includes a high precision tracking system, which 
consists of a silicon vertex detector and dedicated tracking planes. The tracking planes 
consist of silicon microstrip detectors in the region with high charged-particle flux close 
to the beam pipe and straw tube detectors which provide coverage up to the edge of the 
LHCb geometrical acceptance. The tracking planes are located either side of the dipole 
magnet to allow the measurement of the momenta of charged particles as they traverse the 
detector. Excellent particle identification capabilities are provided by two ring imaging 
Cherenkov detectors which allow charged pions, kaons and protons to be distinguished 
from each other in the momentum range 2-100 GeV/c. The energy of particles traversing 
the detector is measured using a calorimeter system which is sensitive to photons and 
electrons, as well as hadrons. Muons are identified using a dedicated detector system. 

The experiment employs a multi-level trigger comprised of a hardware trigger which 
uses information from the calorimeter and muon system and a software trigger which 
performs a full reconstruction of the event, including tracks and vertices. 

The simulated events used in this analysis are produced using the Pythia 6.4 genera- 
tor (13 , with a choice of parameters specifically configured for LHCb [14] . The EvtGen 



package 15 describes the decay of 6-hadrons and the Geant4 toolkit 16 simulates the 



detector response, implemented as described in Ref. 17 . QED radiative corrections to 



the Bg—> K + K decay are generated with the Photos package 18 



sectionTrigger and event selection At LHCb, 6-hadrons are produced with an average 
momentum of around 100 GeV/c and have decay vertices displaced from the primary 
interaction vertex. Combinatorial background candidates, produced by the random com- 
bination of tracks, tend to have low momentum and originate from a primary pp collision 
vertex. These features are typically exploited to select 6-hadrons and reject background. 
The distance of closest approach (impact parameter) of 6-hadron decay products to any 
primary vertex is a particularly important discriminant in the trigger because it is an 
order of magnitude faster to compute than the momenta of the same decay products. 
For this reason, the majority of triggers for hadronic 6-hadron decays begin by selecting 
tracks with a significant displacement from any primary vertex. However, such require- 
ments introduce a time-dependent acceptance which biases the decay time distribution of 
the selected 6-hadron candidates and a significant investment of effort is often required to 
correct for this bias. 

The analysis presented here uses an approach that selects 6-hadrons without bias- 
ing the decay time distribution, other than trivially through a simple minimum decay 
time requirement, limiting the systematic uncertainties associated with correcting for any 
time-dependent acceptance effects. This is achieved using neural networks based on the 



NeuroBayes package 19 in the software trigger and event selection. Neural networks 



have advantages over traditional "cut-based" approaches since they are able to exploit 



2 



the correlations between variables in order to increase signal purity, allowing 6-hadrons 
to be selected without resorting to requirements on impact parameters or flight distance. 

The LHCb software trigger has two stages which run sequentially. Due to restrictions 
on processing time it is not possible to employ a neural network in the first level of the 
software trigger. Instead, only tracks that are not used in the first level decision are passed 
to the second trigger level in order to avoid a potential bias. These tracks are required to 
pass a loose pre-selection with requirements on their momenta, transverse momenta and 
track fit quality. The tracks are then combined to form B meson candidates, using a kaon 
mass hypothesis for both tracks, and further requirements are made on the distance of 
closest approach of the two tracks to each other, the mass of the resulting candidate, the 
helicity angle of the tracks in the B meson rest frame and the quality of the decay vertex 
fit. 

After this pre-selection the candidates pass through a first neural network, trained on 
simulated B° s — > K + K~ , B° — > K + tt~ and background events, which uses the momenta 
and transverse momenta of the tracks and B meson candidate, the distance of closest 
approach of the two tracks, helicity angle, the x 2 °f the vertex fit and the uncertainty 
on the fitted B meson mass to discriminate between signal and background. After this 
stage the data rate is reduced to a level such that each event may be fully reconstructed, 
including information from the particle identification system. A second network, trained 
on the same simulated events, uses the information presented to the first network along 
with particle identification information to further increase the purity of B mesons in the 
selected candidates. 

Roughly half way through 2011 the luminosity delivered by the LHC accelerator in- 
creased to a level such that it was necessary to require that the decay time of B meson 
candidates exceeded 0.3 ps in order to keep the trigger rate within acceptable limits. This 
requirement only biases the decay time distribution in a trivial way, except through a 
possible difference in the decay time resolutions of the trigger and offline reconstruction 
software. 

After the trigger, the tracks associated to the selected candidates are removed from 
the primary vertex fit to avoid a potential bias in the measured decay time. The purity of 
signal candidates is then further enhanced using two additional sequential neural networks. 
The first network is trained using simulated events and combines the same information 
used by the trigger networks along with particle identification information, the energy 
of each track from the calorimeter, the probability that either track is formed from the 
association of random hits in the detector and the x 2 P er degree of freedom for both track 
fits. This network benefits from the more detailed full event reconstruction which is not 
available in the trigger. 

The second network is trained on the data recorded in 2011 using s Weights [20] , 
which are calculated in a window around the signal peak and in the upper sideband 
region (5.45 < itlk+k- < 5.85GeV/c 2 ) of the invariant mass spectrum. The sW eights 
are obtained from a fit to the invariant mass spectrum of the candidates and the neural 
network uses them to discriminate between signal and background. This network uses 
the output of the first network as input, all the input variables used by the first network, 
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the uncertainty on the decay time of the B meson candidate and the impact parameter of 
the B meson candidate with respect to the primary interaction vertex. Only candidates 
with a decay time of r > 0.3 ps are used in the network training. 

The event selection is determined by making a requirement on the output of this second 
neural network that maximises the metric s / a/s + b, where s is the number of signal 
decays in the region 5.05 < m K + K - < 5.85GeV/c 2 and b is the number of background 
combinations. 

The trigger and offline software reconstruct B meson decay times with different reso- 
lutions. Potential "edge-effects" introduced by the trigger requirement that r > 0.3 ps are 
avoided by requiring that candidates satisfy r > 0.5 ps in the final event selection. The 
contribution from the B° — >■ K + n~ and B®—)- K + K~ modes are separated by demanding 
tight requirements on the particle identification properties of the final state particles. A 
small level of contamination from decays of baryons is further suppressed by demanding 
that the final state particles are not compatible with the proton hypothesis. 



3 Analysis of the effective B Q s — > K + K lifetime 

The effective B° s — > K + K~ lifetime is evaluated using an unbinned log-likelihood fit. A 
fit to the invariant mass spectrum is performed to determine the sW eights that are used 
to isolate the B° s — > K + K~ decay time distribution from the residual background. The 
B° s — > K + K~ signal component is described by a Gaussian function. The background 
contamination from partially reconstructed B meson decays is described by a further 
Gaussian function and the combinatorial background is described by a Chebychev poly- 
nomial with one free parameter. It should be noted that the kaon mass is assigned to both 
final state particles in the vertex fit and hence the reconstructed B° — > K + tt~ mass is 
shifted towards higher values than the nominal mass, creating an asymmetric distribution. 
The B° K + ir~ signal component is therefore described by a Crystal Ball function 21 



with the tail on high mass side. The parameters of this distribution are fixed using a fit to 
the independent B° — > K + ir~ sample, separated using particle identification information. 

The fit finds 997 ± 34 B° K + K~ decays and 78 ± 17 B° K + n~ decays in the 
data with 253 ± 25 and 169 ± 20 combinatorial background and partially reconstructed 
combinations respectively. Figure [ija) shows the resulting invariant mass spectrum for 
B®—> K + K~ candidates. 

Using the sWeights returned by the mass fit, the B° s — > K + K~ decay time distribution 



is extracted from data using the sPlot technique 20 . Since there is no acceptance bias 
to correct for, the lifetime is determined using a fit of the convolution of an exponential 
and Gaussian function to account for the resolution of the detector. The mean of the 
Gaussian function is fixed to zero and its width is fixed to the expected resolution from 
simulated events, which is a t = 0.04 ps. 

The effective B® — > K + K~ lifetime is found to be 

t kk = 1.455 ± 0.046 (stat.) ps. 
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Figure 1: (a) Invariant mass spectrum for all selected B® — > K + K~ candidates, (b) 
Decay time distribution of B° s — > K + K~ signal extracted using sWeights and the fitted 
exponential function. 



Figure [l[b) shows the corresponding fit to the decay time distribution of B® — > K + K~ 
signal. 

Since the decay B° — > K + ir~ has similar kinematics, it can be used as a control mode. 
However, since the kaon mass hypothesis is assigned to both tracks, the measured decay 
time is biased to larger values for B° — > K + tt~ . To avoid this bias a fit is made to the 
reduced decay time, which is defined as the decay time divided by the invariant mass. This 
quantity is independent of the mass assigned to the two tracks and is also unbiased by the 
selection, following an exponential distribution with decay constant equal to m B o/r B o. 

Using the value of the B° mass (7j as input, the B° lifetime is found to be 

r B o = 1.536 ±0.031 (stat.) ps 
which agrees with the current world-average t b q = 1.519 ±0.007 ps @. 



4 Evaluation of systematic uncertainties 

A wide range of effects that can influence the measurement of the effective B® — > K + K~ 
lifetime has been evaluated. The individual contributions to the systematic uncertainties 
are described below and their estimated values are summarised in Table [Q 

The key principle of this analysis is that the trigger and event selection do not bias 
the decay time distribution of the selected B® — > K + K~ candidates other than in a 
trivial way through a minimum decay time requirement. This has been tested extensively 
using simulated events at each stage of the selection process to demonstrate that no step 
introduces a time-dependent acceptance. Figure [2] shows the efficiency of the full trigger 
and event selection as a function of decay time for simulated B° s — > K + K~ candidates. The 
graph is fitted with a first order polynomial with a gradient of — 0.09±0.30 ns _1 consistent 
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with a uniform acceptance. Possible discrepancies between simulated and real events are 
considered by comparing the distributions of variables used by the neural networks and 
good agreement is observed. The available quantity of simulated events limits any non- 
zero gradient in the acceptance to within 0.30 ns -1 . This limit is used to evaluate the 
shift in the measured effective lifetime due to the presence of a linear acceptance and a 
negligible deviation is observed and is not considered any further. 




6 8 10 
Decay time [ps] 



Figure 2: Combined efficiency of LHCb trigger, selection neural networks and particle 
identification requirements as a function of decay time for simulated B® — > K + K~ signal 
candidates. 



Studies using simulated events have demonstrated that the efficiency with which tracks 
are reconstructed decreases as the impact parameter of the track with respect to the 
beam line (IP 2 ) increases. This introduces a decay time acceptance that may bias the 
measured lifetime. Such a systematic bias has been evaluated using a combination of 
data and simulated events. First, the effective lifetime of simulated — > K + K~ signal 
candidates is found after reconstruction to deviate by 5 fs from the generated value. 
Second, the tracking efficiency is parametrised as a function of IP 2 using simulated events. 
The calculated efficiency is then applied as a weight to events in data according to their 
IP 2 values and the effective lifetime is evaluated. This produces a deviation of 4 fs with 
respect to the unweighted events. The larger of these two shifts is taken as the systematic 
uncertainty introduced by the reconstruction acceptance. 

The invariant mass distribution of B Q S — > K + K~ signal candidates is modelled using 
a Gaussian function. Potential systematic effects due to this parametrisation are evalu- 
ated by using the sum of two Gaussian functions to model additional resolution effects 
and separately a Crystal Ball function 21 to model final state radiation. Additionally 
the background parametrisation is checked by replacing the first order Chebychev poly- 
nomial with an exponential function. All these changes shift the measured lifetime by 
approximately 1 fs which is taken as the systematic uncertainty. 

The decay time distribution is fitted with an exponential function convolved with 
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a Gaussian function to model detector resolution, where the resolution is fixed to the 
value obtained from simulated events. As a cross-check, the fit is performed with the 
resolution parameter allowed to vary and also using a simple exponential function without 
attempting to model detector resolution. No deviation from the default measurement of 
the effective lifetime is observed in either case. 

The effective B° s — > K + K~ lifetime measurement has been evaluated using an al- 
ternative method which makes a simultaneous fit to the invariant mass and decay time 
distributions. This approach requires a parametrisation of the background decay time dis- 
tribution since the sPlot technique is not used. Both methods give equivalent numerical 
results. 

A wide range of different approaches to the training of the neural network have been 
tested, as well as the influence of different alignment and calibration settings and the 
number of simultaneous primary interactions in the detector. All results obtained in 
these checks are consistent with the result of the default analysis. 

The measured decay times of B meson candidates are determined from the distance 
between the primary interaction and the secondary decay vertex in the silicon vertex 
detector. A systematic bias may therefore be introduced due to uncertainty on the LHCb 
length scale. This effect is estimated by considering the uncertainty on the length scale 
from the mechanical survey, thermal expansion and the current alignment precision. The 
uncertainty on the length of the detector along the beam-line is determined to be the 
dominant effect and a corresponding systematic uncertainty is assigned. 

The effective lifetime is obtained by fitting a single exponential function to the dis- 
tribution given by Eq. [TJ However, the requirement that the decay time be greater than 
0.5 ps diminishes the IY component relative to the Th component in the decay time dis- 
tribution. This effect has been evaluated using simulated events and a deviation of 1 fs 
from the result of a fit to the full decay time range is observed. 

If the production rates, R, of B° s and B® mesons are not equal then an additional 
oscillatory term is introduced into the decay time distribution given in Eq. [TJ proportional 
to the production asymmetry A P = [R(B° S ] - R(B°)] / [R(B° 8 ) + R(B°)] . This term may 
alter the measured effective lifetime. Since the B® meson shares no valence quarks with 
the proton Ap(B° s ) at LHCb is expected to be small. Making the conservative assumption 



that the \A P (B°)\ = \A P (B°)\ = 0.01 [22) we find a shift from the expected value of the 
effective lifetime of 2 fs using simulated events. This value is assigned as the systematic 
uncertainty. 



5 Conclusions 

Two-body charmless B decays offer a rich phenomenology to explore the phase structure 
of the CKM matrix and to search for manifestations of physics beyond the SM. The 
effective lifetime of the decay B® — > K + K~ is of considerable theoretical interest as it is 
sensitive to new particles entering at loop level. A measurement of this quantity is made 
possible by the excellent particle identification capabilities of the LHCb experiment. 
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Table 1: Contributions to the systematic uncertainty on the effective B° s — > K + K~ lifetime 
measurement. The total uncertainty is calculated by adding the individual contributions 
in quadrature. 



Systematic sources 


Uncertainty on tkk [fs] 


Reconstruction efficiency 


5 


Signal model 


1 


Background model 


1 


Length scale 


1 


Minimum decay time requirement 


1 


Production asymmetry 


2 


Total 


6 



The effective lifetime of the decay mode B s — > K + K~ is measured using 1.0 fb~ of 
data recorded by the LHCb detector in 2011. A key element of this analysis is that the 
trigger and event selection selects B mesons without biasing the decay time distribution. 
This is achieved using a series of neural networks. Although this dedicated trigger has a 
lower efficiency compared to the one used in the previous LHCb measurement p], it has 
the advantage of avoiding systematic uncertainties related to the depletion of candidates 
at low decay times and provides an independent approach to measuring the B® — > K + K~ 
effective lifetime. It is measured as 

T KK = 1-455 ± 0.046 (stat.) ± 0.006 (syst.) ps, 

in good agreement with the SM prediction of 1.40 ±0.02 ps and with the measurement on 
data recorded by LHCb in 2010 of 1.440±0.096 (stat.)±0.008 (syst.)±0.003 (mod.) ps |6). 
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